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ABSTRACT 

Context. Multi-band observations of planetary transits using the telescope defocus technique may yield high-quality light curves 
suitable for refining the physical properties of exoplanets even with small or medium size telescopes. Such observations can be used 
to construct a broad-band transmission spectrum of transiting planets and search for the presence of strong absorbers. 
Aims. We have thoroughly characterised the orbital ephemeris and physical properties of the transiting planet and host star in the 
WASP-23b system, constructed a broad-band transmission spectrum of WASP-23 b and performed a comparative analysis with 
theoretical models of hot Jupiters. 

Methods. We observed a complete transit of WASP-23 b in seven passbands simultaneously, using the GROND instrument on the 
MPG/ESO 2.2 m telescope at La Silla Observatory and telescope defocussing. The optical data were taken in the Sloan g', r', i' 
and z' passbands. The resulting light curves are of high quality, with a root-mean-square scatter of the residual as low as 330 parts 
per million (ppm) in the z'-band, with a cadence of 90 s. Near-infrared data were obtained in the JHK passbands. We performed a 
MCMC analysis of our photometry plus existing radial velocity data to refine measurements of the ephemeris and physical properties 
of the WASP-23 system. We constructed a broad-band transmission spectrum of WASP-23 b and compared it with a theoretical 
transmission spectrum of a Hot Jupiter. 

Results. We measured the central transit time with a precision ~8 s. From this and earlier observations we obtain an orbital period 
of P = 2.9444300 ± 0.0000011 d. Our analysis also yielded a larger radius and mass for the planet (R p = 1.067+^^ R Jup and 
M p = 0.917+j555*° Mj„ p ) compared to previous estimates (R p = 0.962+55°56 Rj» P and M p = 0.884+°, 55^4 Mj^). The derived transmission 
spectrum is marginally flat, which is not surprising given the limited precision of the measurements for the planetary radius and the 
poor spectral resolution of the data. 

Key words. Stars: general - Stars: planetary systems - Planets and satellites: general - Planets and satellites: fundamental 
parameters - Planets and satellites: atmospheres 
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observed radius of a planet with a surface gravity of 15 ms~ 2 , 
over the wavelength interval 350-700 nm. For the purposes of 
transmission spectroscopy it is practical to convert these planet 
variations to the change of a directly observed quantity such as 
the transit depth 5 - (R p /R*) 2 . Assuming a radius variation with 
an amplitude « = 0.03 (3%), the transit depth variation is defined 
as 

n(R p + nR v f ttR 2 

AS = p - p; 1 * 2n6 (1) 

nRl nRl 

A transiting Jupiter orbiting a Sun-like star causes 6 = 0.01 
which would translate into a variation of A<5 = 6 x 10 (or 
A(Rp/R„) = y/2nS ~ 0.024). Such precisions may be achieved 
by modelling light curves obtained using ground-based 2-4 m 
telescopes. The best results can be obtained by defocusing the 
telescope, which averages out sensitivity variations between the 
pixels on the detector and improves the observing efficiency by 
enabling a larger fraction of open-shutter time. However, the 
bandwidths of the filters used (which give the spectral resolu- 
tion) is the other critical factor affecting the chances of detecting 
the signal in Equation 1, similar to the detection of spectral fea- 
tures in spectrophotometry. 

Several successful searches for the theoretically predicted 
strong absorbers in the optical and near-infrared (NIR) have been 
conducted first from space using the Hubble and Spitzer space 
telescopes, resulting in the dete ction of Na, CO, H?Q, CH 4 and 
Rayleigh scattering in the blue (Charbonneau et al. 2002; Pont 
et al. 12008b ISing et al.ll2011bh . Hampered bv the atmosphere 
of the Earth (e.g. telluric contamination) and instrumental sys- 
tematics, the ground-based detection of spectral features in ex- 
oplanet atmosphere s took longer to yield re sults (iRedfield et al.1 
20081 ISnellen"etani2008t ISing et aUGOl iaLl2012h . 

Although transit photometry or spectrophotometry is often of 
high quality, the method pushes the limits of current instrumen- 
tation and has yield ed controversial results (ISwain et al.l 12008: 
iBeaulieu et al.ll20Tlb . The success of the method may also be 
hampered by the activity of the host star, in particular spot cross- 
ing events which can significantly increase the complexity of 
light curve analysis and lower the prec ision of the planetary ra- 
dius measurement dKnutson et alj|20l"ll) . 

The transiting hot Jupiter WASP-23 b orbits a moderately 
bright (V = 12.68) Kl V (T eff = 5150+ 100 K) s tar and was 
found bv the Super WASP survey (Tria udet al. 1201 ft Pollacco et 
al. |2006J). The planet (radius 0.96 + 0.05 R h mass 0.88 ± 0. 1 Mj) 
has a 2.94 d circular orbit (e < 0.062 at the 3<x level). 

In this paper we present ground-based simultaneous optical 
and NIR photometry of a transit of WASP-23 b. We refine its 
measured properties and orbital ephemeris. We then construct 
the first transmission spectrum of WASP-23 b and investigate 
the variation in radius of the planet with wavelength. The lay- 
out of the paper is as follows. In Section 2 we present the ob- 
servations. Section 3 discusses our data reduction methods. In 
Section 4 we report the analysis of the data, measure the physi- 
cal properties of the system, construct the planet's transmission 
spectrum, and compare this theoretical predictions. We conclude 
in Section 5. 



2. Observations 

We monitored the flux of WASP-23 (also known as GSC 07635- 
01376, a = 06 h 44 m 30.65 s , 6 = -42°45'41.0"; J2000.0; V = 
12.68 mag) during one transit on February 25, 2012 (Fig. 1 & 
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Fig. 1. Simultaneous optical to near-infrared (Sloan g',r',i','z and 
Johnson J, H, K) raw transit light curves of WASP-23 (dots), obtained 
with GROND on UT February 25, 2012 an the corresponding correction 
functions (red lines). 

2). The data were collected using the GROND instrument, at- 
tached to the MPG/ESO 2.2m telescope at ESO La Silla Ob- 
servatory (ChilefJ. GROND (Gamma Ray Burst Optical and 
Near-Infrared Detector) is a simultaneous multi-channel imager 
that has been specificall y designed for gamm a-ray burst (GRB) 
afterglow observations (IGreiner et al. 1120081) . Its capability to 
simultaneously monitor sources in seven optical to NIR pass- 
bands make it an excellent instrument for obtaining high qual- 
ity, multi-band transit light curves. It is capable of yielding 
light curves with a root-mean-square (r.m.s.) scatter better than 
1 mmag, depending on the atmospheric conditi ons dLendl et al.l 
l201(tlNikolov et al. 11201 2l : lMancini et alJl2013h . 

To reach GROND's detectors the light from the telescope 
is split using dichroics into seven beams, resulting in response 
curves close to the Sloan g' (A = 459 nm), r' (A = 622 nm), 
i' (A = 764 nm), z' (A = 899 nm) in the optical and J (A = 
1256 nm), H (A = 1647 nm) and K (A =2151 nm) bands in 
the NIR, respectively. The optical beams are incident onto four 
2048 x 2048 back-illuminated E2V CCDs with 5.4' x 5.4' fields 
of view (FOVs) at a plate scale of 0".158 pixel" 1 . The NIR 
channels are imaged onto three 1024 x 1024 Rockwell HAWAII- 
1 arrays with 10 x 10 arcmin FOVs and plate scales 0".6 pixel" 1 . 
The off-axis guide camera was used to autoguide the telescope. 

1 Detailed information on the 2.2m telescope, the GROND in- 
strument and the relevant technical settings can be found at 
http://www.eso.org/sci/facilities/lasilla/telescopes/2p2/ 
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Fig. 2. Le/f panel: Detrended transit light curves (dots) of and the best-fit models (lines) for the WASP-23 GROND data. The optical photometry 
clearly exhibits a change in the shape of the transit due to the different limb-darkening in each band. Right panel: Residuals of the light curve 
compared to the fitted models. 



Light curves were acquired in GROND 's non-dither mode to 
decrease systematics associated with inter-pixel sensitivity vari- 
ations. We obtained data for WASP-23 and eleven and seven 
reference stars in the optical and NIR images, respectively. We 
applied a moderate defocus to the 2.2m telescope, such that the 
point spread functions (PSFs) had diameters ~ 76 pixels (12") in 
the optical and ~ 13 pixels (7. "6) in NIR. Table 1 gives details 
of the instrumental setup and illumination levels of the seven 
beams. As the guide camera was also defocussed, we made sure 
in advance that at least one bright star was within its FOV. 

The main advantages of the telescope defocu s technique 
have been discussed by s e veral researchers, e.g. [ Winn et al.l 
(l2007h . lGillon et all (l2007h ISouthworfh et al. I (120101) (and refer- 
ences therein). ICroll et al71 (l2010h . Although the method results 
in non-uniformly defocused PSFs of the stars across large FO Vs 
(> 1° x 1°), the technique performs well on small FOVs (such 
as GROND). It allows superb PSF sampling, significantly reduc- 
ing the random and systematic errors resulting from detector im- 
perfections and resposne variations. The technique also avoids 
detector saturation due to fluctuations in seeing, because such 
fluctuations are small compared to the level of defocusing. We 
note that telescope defocusing is not well suited to photometry 
in crowded fields or when the target star has nearby companions. 
We used the 2MASS catalogue to rule out these possibilities for 
our observations of WASP-23. 

We observed WASP-23 on February 25 2012 for 4.4 hr from 
UT 00h32m to 04h56m. An exposure time of 90 s was used in 



Table 1. Instrumental and observational settings used to obtain data for 
WASP-23 (median values), during the observation on UT February 25, 
2012. 



GROND 
band 


Gain 
(ADU) 


LI nn g 

(pixel) 


Max count 
(ADU)xlO 3 


Background 13 
(ADU)xlO 3 


g' 


1.45 


57 


18.6 


0.3 


r' 


1.33 


67 


43.8 


0.3 


r 


1.62 


51 


27.2 


0.3 


z' 


1.74 


48 


22.3 


0.4 


I 


2.39 


9 


26.7 


5.0 


H 


2.65 


6 


18.2 


17.0 


K 


2.55 


24 


13.1 


11.0 



Notes. (a) Diameter of the defocussed psf ring in pixels. (b) Measured 
using an annulus, centered on the stellar psf, encompassing an area four 
times larger than the area used to produce a photometric measurement 
of the stellar psf. 



the optical and resulted in 91 integrations. The slow read-out 
mode was selected (~ 48 s read-out- time) as it delivers high- 
quality images Q. In the NIR averaged stacks of seven 10 s 
integrations were obtained to improve the signal-to-noise ratio 

2 The GROND detector can also be read-out in 'fast' mode (~ lis 
read-out-time), but it has been found during previous transit observa- 
tions that images obtained in that mode suffer from systematic effects 
that lower the quality of the resulting light curves. 
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Table 2. r'-band light curve extract 



BJD a 


Flux 


0"Flux 


z c 


x d 


y e 


2.52787 


1.000410 


0.000256 


1.03539 


7.73 


71.17 


2.52990 


1.001160 


0.000269 


1.03370 


7.61 


70.30 


2.53193 


0.999446 


0.000253 


1.03274 


7.54 


70.06 


2.53396 


1.000615 


0.000255 


1.03170 


6.95 


70.66 


2.53599 


0.998938 


0.000253 


1.03100 


6.72 


70.67 


2.53802 


1.000043 


0.000252 


1.03079 


6.18 


69.93 



Notes. (a) BJD - 2455980; (b) 1-cr error-bar of the flux, as measured by 
APER; (c) airmass; (d) x - 630 horizontal pixel position of the defocussed 
psf; <d) y - 1300 vertical pixel position of the defocussed psf; 



(SNR), resulting in a total of 181 images. At the beginning of 
the observation the airmass of the target began at 1 .04, peaked at 
1 .03, and increased to 1 .47 by the end of the observing sequence. 
The observing conditions were nearly photometric, with seeing 
variation from ~ 16 to ~ 18 pixels as measured from a fit to the 
wings of the defocused PSFs. 

Theoretical models of exoplanet transmission spectra imply 
that the NIR region is richer than the optical in terms of po- 
tential spectroscopic signatures. We therefore aimed to obtain 
high-precision light curves for that region as well. The control 
of GROND allows only a single integration time (limited from 
to 10 s) for the three JHK bands simultaneously, so we opti- 
mised the exposure time such that the highest SNR was obtained 
for the band when the star had the most flux (J). For the degree 
of telescope defocus, we aimed to optimise the SNR in the op- 
tical, as we expected better performance in that region based o n 
previous experience (Ni kolov et al. II20T2I; iMancini et al.ll2013h . 

3. Data reduction and analysis 

The data were reduced and analysed on a channel-by-channelba- 
sis, using a customised pipeline. GROND's optical and NIR data 
are provided in separate fits files. The images from the optical 
channels are packed into four-extension data cubes and the JHK 
read-outs are delivered as a single 3072 x 1024 image. After pro- 
cessing to collapse the optical data cubes and to split each NIR 
image, the data were sorted into science and calibration images. 

For the optical data a median-combined master bias frame 
was computed using 20 zero-second images. We estimated 
the average dark current of GROND's optical detectors to be 
~ O.Ole s 'pixeF , which is negligible for our data. We there- 
fore neglected the dark current correction. A median-combined 
master flat-field was calculated using six bias-corrected and 
dithered twilight sky flat-fields that survived a selection crite- 
rion for linearity in each of the passbands. Finally, each optical 
science frame was de-biased and flat-field corrected to produce 
calibrated time series of images. 

In the NIR, a series of ten dark frames with integration time 
of 10 s were employed to compute a median-combined master 
dark frame. A series of 24 dithered twilight sky flat frames were 
used to obtain a median-combined master sky flat. The mas- 
ter dark frame was subtracted from all the science and sky flat 
frames. Each dark-subtracted science image was then corrected 
with the master sky flat frame. An electronic odd-even readout 
pattern along the y-axis was removed before flat-fielding. This 
was done by smoothing each image and comparing it with the 
unsmoothed one after the master dark had been subtracted. The 
amplitudes of the read-out pattern were then determined by com- 
paring the median level of each column to the overall median 



Table 3. Theoretical limb-darkening coefficients, light curve residual 
scatter (in parts per million, ppm) and the /J-factor used to estimate the 
photometric uncertainties of the GROND data. 



c Dnwn 
(jKUIN d 


U 1 


U2 


r.m.s 


Q 

P 


band 






(ppm) 




g' 


0.7413 


0.0719 


419 


1.00 


r' 


0.5323 


0.1804 


371 


1.00 


i' 


0.4224 


0.1983 


430 


1.27 


z' 


0.3491 


0.2062 


330 


1.28 


J 


0.1063 


0.2589 


1480 


1.50 


H 


0.1126 


0.3130 


1760 


1.52 


K 


0.2389 


0.2467 


2980 


1.16 



level. Each column was then shifted to the overall median level. 
The master sky flats were divided out from the science images 
after removal of the read-out pattern. 

Previous investigators have commonly obtained dithered im- 
ages prior to the main observing sequence to further correct the 
science images against sky background variations. However, 
similar corrections are inadequate for sky gradients that vary on 
timescales less than the duration of the observations. In addition, 
relative photometry, if performed at gradient-stable sky back- 
grounds is in practice insensitive to sky variations that concern 
simply the sky level. Finally, previous analyses of WASP-44b 
(Man cini et al.ll20l3h with and without sky corrections resulted 
in indistinguishable light curves with similar scatter. We there- 
fore did not obtain dithered sky images. 

We performed circular aperture photometry on each of the 
calibrated images. To accurately determine the central positions 
of the heavily defocused PSFs, we proceeded as follows. First, 
we extracted sub-images, roughly centred on each star. The 
brightest stars were rings with typical diameters ~ 50 and ~ 15 
pixels in the optical and NIR, respectively. We then convolved 
each image with a Gaussian kernel to produce a smoothed PSF 
with a well-defined peak, located at the centre of the original 
PSF. Finally, we fitted a two-dimensional Gaussian to measure 
the PSF centre and full width at half maximum, which was lat- 
ter used to monitor the seeing variations during the observation. 
Instrumental fluxes and uncertainties (computed taking into ac- 
count a random photon noise) were acquired usinga daophot- 
type photometry algorithm as implemented in the idi0 apeeQ rou- 
tine. 

To compute the relative flux light curve of WASP-23 in 
the four optical channels we used 2MASS J06444546-4242083 
as a comparison star. Because of the relatively small FOV of 
GROND we performed an exhaustive investigation of the possi- 
ble reference stars prior to the observations. This object satisfied 
the requirements, being isolated enough to avoid contamination 
with bright neighbours when defocused, and being the only stel- 
lar source in the optical FOV of a similar brightness and colour 
to WASP-23, thus reducing systematic effects due to differential 
colour extinction. For instance, the instrumental magnitude dif- 
ferences (WASP-23 minus comparison star), measured on UT 
Februray 25 2012 were Ag' = 0.5 10 mag, Ar' = 0.082 mag, 
A/' = 0.3 19 mag, Az' = 0.785 mag. For the same reasons, in 
the NIR channels we used 2MASS J06445925-4245061 as the 
comparison star. The instrumental magnitude differences were 
AJ = -0.784 mag, AH = -0.696 mag and AK = -0.396 mag. 

3 The acronym IDL stands for Interactive Data Language. For further 
details see http://www.exelisvis.com/ ProductsServices/IDL.aspx 

4 part of aper Astronomy User's Library. For detailed information see 
http://idlastro.gsfc.nasa.gov/ 
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To select optimum sizes for the inner aperture and sky an- 
nulus we investigated the root-mean-square (rms) of the out-of 
transit baseline flux as a function of various combinations of 
aperture and annulus sizes. Table 2 exhibits an r'-band light 
curve extract along with the airmass and measurements of the 
(x, y) drift of the PSF on the detector employed to correct the 
light curve during the analysis. 

We used the UT time stamps available in the fits headers of 
the data and computed the central times of each integration using 
the exposure times. All times were then converted to Barycentric 
Julian Date (BJD), using the converter of Eastman et al. (2010). 
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Fig. 3. Observed minus computed (O-C) transit times for WASP-23 b, 
using observations reported in the literature (black dots) and our mea- 
surement (orange dot), compared to a linear ephemeris (dashed line). 



4. Light curve analysis 

We performed the light curve analysis in two steps. First, we 
simultaneously fitted the seven optical to NIR GROND light 
curves to derive accurate estimates of the transit central time T c 
and orbital period P, and evaluated the orbital parameters and 
planetary radius (/, a/R*, R p /R*). 

In the second step we complemented the GROND photom- 
etry with t he radial velocity (R V) measurements of WASP-23, 
reported bv lTriaud et al. I (I201 ll) and perfor med a joint photome - 
try plus RV MCMC analysis using exofast dEastman et alJ2 012). 
As our ultimate goal is to construct a ground-based transmission 
spectrum of WASP-23, based on uniform transit parameters from 
all bands we excluded the RVs obtained during transit, which are 
affected by the Rossiter-McLaughlin effect. 

As exofast allows an analysis of one light curve at a time 
we use the results from the first step (the simultaneous fit to the 
photometry) as priors for exofast. The latter was used to derive 
reliable estimates of the physical properties and their uncertain- 
ties for WASP-23 b a nd its parent star, e xofast utilizes empirical 
polynomial relations dTorres et al.l2 010) between the masses and 
radii of stars, and their surface gravities (log g), effective temper- 
atures (T e s) and metallicities [Fe/H], based on a large sample of 
well-studied non-interacting binaries. Including this empirical 
relation as a constraint allows a derivation of the full orbital and 
physical properties from transit light curves and RV data (East- 
man et al. l2012l) . 

In the first step, we modelled the transit light curves simul- 
taneously in flux units). A two-component correction function 
was constructed, containing a transit model multiplied by a base- 
line function that accounts for the most significant known sys- 



Table 4. Derived system parameters for WASP-23 b from light curve 
fits, compared to the corresponding results reported in the discovery 
paper. The final estimates of our analysis are the weighted mean values 
from all bands. 



Band 



GROND g' 
GROND r> 
GROND V 
GROND zf 
GROND J 
GROND H 
GROND K 



R P /R. 
013402™ 52 
0.13336 



1343 3 +0 
u.ijtj -ao 0043 

13400+ 000042 
1340 +(Vo( ' 12 

13520+ 00012 
u.ijj^u p 2 

0.1326 +0 * 14 



Weighted mean 



87.93+J^ 

87.85!| 
88-02+1 

87.83+| 
88.12+0* 

87.75!g 
88.74 H 



a/R* 
9 900+ 014 
9.850+ 014 



9.910 
9.850 



Triaudetal. (2011) 



1 3iQ4+W-(I)()52 
"• iJ " w -n.nnn45 



87.91 



m 

-0.16 
+0.15 

I Q 74+8:4+ 
10.59+i| : |[ 



9.893 



88.39!»-£ 



Tfffffr 

0070 



10.61 



tematic effects affecting ground-based transit observations: air- 
mass and detector intra-pixel sensitivity variations (Fig. 1). 

The first component is based on the analytical transit models 
of iMandel & Agoll d2002l) . which in addition to the transit cen- 
tral time (T c ) and orbital period (P) is a function of the orbital 
inclination (/), planet semi-major axis and radius normalized to 
the stellar radius (a/R*) and (R p /R t ), respectively. To account 
for the limb darkening of the star we adopted the quadratic limb 
darkening law: 



^ = 1 



Ml(l - jj) - M 2 (l -/i) 



(2) 



where / is the intensity and fi is the cosine of the angle be- 
tween the line of sight and the normal to the stellar surface. We 
fixed the two limb darkening coefficients (u\ and u%) to their the- 
oretical values (see Table 3 for details), which we obtained by 
inter polating within the calculated and tabulated A TLAS mod- 
els dClaret & Bloemenll20lU lEastman et al.ll2012r) to the stel- 
lar parameters T eS = 5150 + 100 K, log g = 4.4 ± 0.2 (cgs), 
[Fe/H] = - 05 + 0.13 and v t = 0.8 + 0.3 km s" 1 , which we 
adopted from lTriaud et al. I (1201 lb . 

We constructed a baseline function to account for systematic 
effects that concern the (in- and out-of-transit) photometry. It 
contains a linear airmass (z) term as well as a quadratic polyno- 
mial of the PSF location (x, y) on the detector to take into account 
for any inter-pixel sensitivity variations: 

F (z, x, y) = 1 + Go + ci\z + ci2X + a^x 2 + a$y + a$y 2 + a^xy. (3) 

To derive the best-fit parameters we employed the 
Levenberg-Marquardt least-squares algorithm, minim ising the 
X 1 function over the data (Fig. 2) (Markwardt 2009). Conve- 
niently, exofast also takes into account the systematics, incorpo- 
rating a baseline function and trends, as described in Eastman et 
al. (120121). We therefore included function (3) in exofast to take 
into account the systematics and to enable a uniform analysis of 
the in- and out- of-transit photometryQ 

An accurate derivation of system parameters from transit 
light curve fits requires a realistic estimation of the related pho- 
tometric errors. The photometric uncertainties produced by 



5 We choose to employ mpfit in the first step as it allows one to con- 
struct a joint model incorporating the entire dataset. In contrast exofast 
permits an analysis on a per-channel basis, i.e. one band at a time. We 
therefore, determined universal values for a/R,, i, T c , P and R p /R, using 
the complete dataset and used these results as priors in exofast, where 
the fitting is done for each band. 
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idl/aper are usually underestimates, so we ignored them in our 
analysis. These are formal statistical uncertainties and do not 
take into account time-correlated noise. We used the out-of- 
transit data of each light curve before detre nding and applied 
the 'time-averaging' procedure proposed by_ Pont et all (|2006) 
and used bv various authors including Gillo n et al.l (120061). Winn 
et al. d2007ll2008h iGibson et alJ (l2008h.lNikolov et al. I (l2012h . 
ISouthworth et all d2012allbl) and lMancini et al.l (1201 2|) to inflate 
the photometric errors to more realistic values. In summary, for 
each band we computed the /?-ratio between the standard devia- 
tion of the out-of-transit unbinned and binned data with bin sizes 
similar to the transit ingress/egress duration. The photometric 
uncertainties were therefore set to the out-of-transit unbinned 
standard deviation multiplied by the corresponding /3-factor (see 
Table 3). 



Table 6. Transit mid-times and O-C residuals computed using tran- 
sit light curves from the literature and the observation reported in this 
work. Note the high precision (~8s) of the the transit time derived using 
GROND's light curves. 



Epoch 


Central time 
(BJD TDB ) 




O-C 

(days) 


Refernce 





2454813.68163" 


00070 


0.000044 


1 


98 


2455102.23664" 


8888 

m 


0.000916 


1 


143 


2455234.73568° 


0.000607 


1 


251 


2455552.73321° 


88818 
88812 

(inns i 


-0.000300 


1 


254 


2455561.56610° 


-0.000700 


1 


397 


2455982.620315°™ 


0.000028 


2 



Notes. 1- Triaud etal. I d201lh : 2- this work. 



4. 1 . System parameters 

To derive system parameters for WASP-23, including the priors 
for exofast, we relied on the light curves from the seven GROND 
passbands. We assumed that i and a should not depend on the 
observed band and constrained them to a single universal value. 
In contrast, we treated the planetary radius as a free parameter 
for each passband. We also set the linear and the quadratic limb 
darkening coefficients to their theoretical values as well as the 
initial guess for the orbital period to the value reported in Triaud 
et al. 42011D. We found i = 88.02° ± 0.19°, a/R, = 9.94 ± 0.11 
and ratios of the radii for each passband: R p /R t (g') = 0.13410 ± 
0.00066, R P IR*</) = 0. 13348 ± 0.00049, R p /R*(?) = 0.13455 + 
0.00047, R p IR»{z') = 0.13439 ± 0.00048, R p /R t (J) = 0.1340 ± 
0.0016, R P /R*(H) = 0.1352 ± 0.0016, R P /R*(K) = 0.1325 ± 
0.0020. These values were input in exofast as priors along with 
the star's 7^, loggr and [Fe/H], fitting the transit photometry 
and RVs simultaneously in each band. The results for the main 
physical properties are summarised in Table 4 & 5. Our result s 
are slightly different from these reported in Triau d et al. 1(1201 ll) . 
though with higher precision, favouring a more inflated planet. 
In particular we improved the precision of i by factors of ~4 and 
~6.5 for the lower and upper limits, respectively. We find the 
radius of WASP-23b to be R p = 1.067+° ^ /Bp , which is larger 
than the radius presented in the discovery paper. 

4.2. Ephemeris 

One of our main goals was to secure the transit ephemeris (T c 
and P) of WASP-23 b. We fixed all physical parameters to their 
values derived in the previous section and fitted GROND's light 
curves simultaneously to derive a transit time. We added this to 
all reported central transit times in the literature (Table 6), then 
performed a linear fit as a function of the observed epoch (E): 



T C (E) = T + E x P. (4) 

We derived P = 2.9444300 ± 0.000001 Id and T c = 
2454813.68159 ± 0.00042 BJD. This peri od is in poor agree- 
ment (~3cr level) with the value reported in Triaud et al. (201 1 ) 
(P = 2.9444256^-°°™°" d), but is more reliable due to the larger 
time span of transit times we used. We investigated the ob- 
served minus calculated (O-C) residuals of the timings versus 
our ephemeris. We found no systematic deviation and conclude 
that there is no compelling evidence for transit timing variations 
(Fig. 3). 



4.3. Radius Variation 

A set of simultaneous multi-band transit light curves provides an 
opportunity to investigate the behaviour of the measured plane- 
tary radius with wavelength. Atmospheric absorption features in 
the optical regime due to sod ium (Na) and potassium (K) have 
been theoretically predicted (Sea ger & Sasselovll2000l) and ob- 



0.138 




Wavelength (|0.m) 



Fig. 4. Radius variation (measured with R,,/R, and indicated with 
yellow disks with error bars) of WASP-23 b in the seven optical to NIR 
passbands, compared to a synthetic transmission spectrum based on a 
planet-wide pressure-temperature profile (grey continuous line). The 
open red diamonds indicate the expected theoretical values, integrated 
over GROND's filter curves. The latter are displayed as dotted lines at 
the origin of the graph. 



0.137 




Wavelength (|xm) 

Fig. 5. Same as Figure 4, but displaying a closer view of the radius 
variation in the optical bands. 
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Table 5. Stellar and plan etary parameters for WASP-23 derived from light curve and RV fits to the GROND photometry and RV measurements 
from iTriaud et al. I (1201 lh . 



Quantity 


Symbol 


This work 


Triaudetal. (2011) 


Stellar parameters 


Mass ( Mq) 


M» 


0.842-;- 


0.78-- 


Radius ( Rq) 


R* 


0.819!; 
Q 421 +0.053 

o i i o+SiStf! 

' 1 -Q.Q44 
4 n 1+0.025 
t.jji _ 0(J29 




Luminosity ( L ) 


1. 




Density (cgs) 
Surface gravity (cgs) 


p, 
\og(g t ) 


1 84^+0.025 
i.otj_ 027 

4.4 + 0.2 


Effective temperature (K) 




5 149 J 


5150 ± 100 


Metallicity 


[Fe/H] 




-0.05 ± 0.13 


Planetary parameters 


Semi-major axis (AU) 
Mass ( Mj) 
Radius (Rj) 
Density (cgs) 
Surface gravity 


a 
M P 
Rp 
Pp 
logigp) 


03798 +000075 
u.uj / yo 00()75 

q 17+0.040 
u - yL -0.039 

1 067 +0045 

-QQ3.8 
0916 +0.036 

-(MB 
o 9 Q4+0.027 

1152+ J 
0770 +() So33 

u.u/ /u A.0035 


0.0376!;^ 
0.884i5f 

962 +0 047 
v.yvz. 0056 


Equilibrium temperature 3 (K) 
Safronov Number 


T 






Incident flux (10 9 erg s _! cm -2 ) 


(F) 


400 + ° M > 




Orbital period (day) 


p 


2.9444300 ±0.0000011 


2 Q444256+ 0000011 


Reference transit time (BJD) 


To 


2454813.68159 + 0.00042 


2455320.12363!]]™; 2 


RV parameters 


RV semi-amplitude (m/s) 


K 


145.57^ 

001039 +oocfo()27 
v-wLvjy _ n .nono27 




Mass ratio 


Mp/M, 





Notes. (a) Assumin g a litt le redistribution over the surface of the planet, i.e., the area that is reemitting with average temperature T eq is 2nR 2 p , as in 
lHansen & BarmarJ l l2007» . 



served in the transmission spectra of hot Jupiters, e.g. XO-2b 
(Sing et al.l 1201 lbl 1201 2|) from the ground. Contrary to this 
example of an agreement between theory and observation, an 
analysis of HST ul traviolet-optica l spectrophotometric data of 
HD 189733 b bv lSing etal.1 d2011ah revealed a featureless spec- 
trum, which might be an indication for the presence of high al- 
titude clouds. These observational results indicate the diversity 
of exoplanet atmospheres, the understanding of which requires 
more observational results. 

To investigate the radius variation of WASP-23 b in 
GROND's passbands we fixed P, i and a/R* to their final results 
from the analysis in Section 4.1, and the linear and quadratic 
limb darkening coefficients to their predicted values. We fitted 
for the planetary radius (specifically the ratio of the radii R p /R t ) 
and the instrumental systematics. This approach improves the 
precision of the measured planetary radius, by avoiding sources 
of uncertainty common to all measurements of R p /R*. In ad- 
dition, the fit takes into account the information for the system 
parameters and reflects differences in the transit depth that are 
determined only by the size of the planet in each passband. To 
estimate realistic values for the normalised planet radii and their 
error bars we carried out the fitting using an MCMC approach. 
The results for the planet radii, obtained at this step are displayed 
in Figs. 4 and 5. The vertical error bars comprise the 68.3% 
confidence levels from the MCMC distributions of R p /R t . As 
expected, the radius measurements in the optical regime show 
much smaller error bars compared to the NIR passbands. This 
is due to the lower scatter of the data in the optical. We also 
show horizontal bandwidths (plotted with error-like bars) that 
have been computed as the full widths at half maximum, using 
GROND's filter curves, which are indicated at the lower part of 
the plot. 



To compare our observational results with theorical predic- 
tions we employed the system parameters from Table 5 and com- 
puted a 1 D model atmosphere of WA SP-23 b, using the code de- 
scribed in Fortnev et al. (20051 120081) . The fully non-gray model 
uses the chemical equilibrium abundances of Lodders & Feg- 
ley (2002) and the opacity database described in lFreedman et al.l 
d2008h . The atmospheric pressure-temperature profile simulates 
planet-wide average conditions. We computed the transmission 
spectrum of the model using the methods described in Fortney et 
al. (l2010l) . For each of the GROND passbands we computed the 
expected theoretical value given the synthetic transmission spec- 
trum, integrating over the corresponding filter curve. Finally, we 
used the derived theoretical values for the planet radius and fitted 
them as a model to the measured values of R p /R t from the seven 
GROND light curves (Fig. 4 & 5). We find a good agreement 
between the observationally derived measurements and synthetic 
values of the planetary radii, consistent with a featureless trans- 
mission spectrum at the ~lcr level. Two of our radius estimates 
(r'-band in the optical and //-band in the NIR) were found to be 
off by slightly more than l<x from their corresponding theoret- 
ical values. This is most likely due to systematics of unknown 
origin. 

5. Summary and discussion 

We have presented observational results based on data collected 
with the GROND instrument on the MPG/ESO 2.2m telescope at 
ESO La Silla, during one transit of the hot Jupiter WASP-23 b, 
the first photometric follow-up of this object after its discovery. 
We employed the telescope defocusing technique and autoguid- 
ing to optimise the quality of the light curve and minimise the in- 
strumental systematics in the complete optical to NIR capability 
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of the GROND instrument for the first time. Light curves con- 
structed from this data set showed a scatter versus fitted models 
as low as 330 ppm (Table 3). For comparison, a recen t observa- 
tion performed bv lTregloan-Reed & Southworthl(l2012l) with the 
3.58 m New Technology Telescope (NTT), in defocussed mode 
obtained light curve scatters of 258 and 211 ppm, during two 
transit events of WASP-50b. Their lower scatter is expected 
due to the larger aperture of the telescope. 

We used the full potential of GROND's optical to NIR 
dataset to refine the ephemeris and system parameters of the 
WASP-2 3b system and c ompa red our results with these pre- 
sented by iTriaud et al. I (1201 lh . We complemented our pho- 
tometric data with existing RV measurements of the parent star 
to derive the full physical properties of the system. We measured 
the central transit time with a precision ~8 s, complemented this 
measurement with published transit epochs, and refined the or- 
bital period of WASP-23 b. We analysed the O-C diagram of the 
transit times and found no significant evidence for transit timing 
variations. Using all of the seven GROND bands we determined 
a larger and more precise value for the radius of WASP-23 b, 
which is in agreement at the ~2cr level with that from Triaud et 
al. (1201 lh . The latter however, was derived based on a joint 
fit to five lower-quality light curves covering the z, R and I+z 
passbands. Despite our larger value for the planetary radius, 
WASP-23 b still sits below the region of the so-called heavily 
"bloated" hot Jupiters, i n accord to th e theor etical models for 
irradiated giant planets of Baraffe et al. (2008). Moreover, Fort- 
ney et al. (2007) presented theoretical mass-radius relationships 
for irradiated giant exoplanets, parametrized by mass, age, and 
an effective orbital distance, defined as the distance from the 
Sun, where a hypothetical planet would receive the same flux 
as the actual planet. Given the luminosity of WASP-23 and the 
semi-major axis of its planet we evaluated that distance to be 
a = 0.585 AU. We then interpolated the models of Fortney et 
al. d2007h to match the distance and the mass of WASP-23b for 
a solar composition at 4.5 Gyr and evaluated R p = 1.042 Rj up , 
which is within the error of our result. 

Finally, we used the potential of the GROND light curves to 
investigate the (possible) variation of the planetary radius with 
wavelength and constructed a ground-based simultaneous opti- 
cal to NIR broad-band transmission spectrum of WASP-23 b. 
According to our results for the system parameters (Table 5) 
the equilibrium temperature of WASP-23 b is T eq = 1 K, 
which places the planet among the pL-class (cooler) of giant 
close-in exoplanet atmospheres in the scheme proposed bv Fort- 
ney et al. (120081) . Taking into account these results, we gener- 
ated a synthetic atmospheric transmission model of an irradiated 
hot Jupiter, assuming a planet-wide pressure-temperature profile 
and compared theory with observations. The optical part of the 
synthetic transmission spectrum at this equilibrium temperature 
and surface gravity (20.73 m s~ 2 ) suggests a radius variation due 
to the presence of gaseous sodium and potassium in the atmo- 
sphere of WASP-23 b, while the NIR parts are dominated by 
water vapour features. The detection of such radius variation 
should be a realistic task, if observations could be taken at the 
precision levels attained in this work and with a higher spectral 
resolution. However, even the narrowest passband of GROND 
(/'-band) spans ~120nm, which effectively reduces the sensitiv- 
ity of our measurements to radius variations. This effect is illus- 
trated in Figs. 4 and 5, where using the filter curves of GROND 
we computed the theoretically expected values for the radius of 
WASP-23 b, given the synthetic spectrum. We fitted our ob- 
servational result for the radius to the theoretical estimates. Al- 
though we found the r'-band, which covers the expected sodium 



absorption line, to deviate from the predicted theoretical value, 
overall the theoretical model is in very good agreement with ob- 
servational transmission spectrum at the ~lcr level. 
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